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Abstract — A fundamental task in wireless communication is chan- 
nel estimation: Compute the channel parameters a signal undergoes 
while traveling from a transmitter to a receiver. In the case of 
delay-Doppler channel, a widely used method is the matched filter 
I algorithm. It uses a pseudo-random sequence of length N, and, 
>— «. in case of non-trivial relative velocity between transmitter and 
-vi receiver, its computational complexity is 0(N 2 logiV). In this paper 
we introduce a novel approach of designing sequences that allow 
feJQfaster channel estimation. Using group representation techniques we 
3 construct sequences, which enable us to introduce a new algorithm, 
called the flag method, that significantly improves the matched 
filter algorithm. The flag method finds the channel parameters in 
C*~) 0(m • N logiV) operations, for channel of sparsity m. We discuss 
O^l applications of the flag method to GPS, radar system, and mobile 
communication as well. 

t-H Index Terms — Channel estimation, time-frequency shift problem, 
h- 1 fast matched filter, flag method, sequence design, Heisenberg-Weil 

^ sequences, fast moving users, high-frequency communication, radar, 

Q GPS. 
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I. Introduction 

A BASIC step in many wireless communication protocols 
lfl6ll is channel estimation: learning the channel parameters 
a signal undergoes while traveling from a transmitter to a receiver. 
In this paper we develop an efficient algorithm for delay-Doppler 
(also called time-frequency) channel estimation. Our algorithm 
provides a striking improvement over current methods in the 
presence of high relative velocity between a transmitter and a 
receiver. The latter scenario occurs in GPS, radar systems, mobile 
communication of fast moving users, and very high frequency 
(GHz) communication. 

Throughout this paper we denote by W = C(Zjv) the vec- 
tor space of complex valued functions on the set of integers 
Zyv = {0, 1, N— 1} equipped with addition and multiplication 
modulo N. We assume that N is an odd prime number. The vector 
space H is endowed with the inner product 



(/i,/a>= E fiMhH 

f° r fiif-2 € H, and referred to as the Hilbert space of (digital) 
sequences. 
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Let us start with example. 

A. Example: The GPS Problem 

A client on the earth surface wants to know his/her geo- 
graphical location. The Global Positioning System (GPS) was 
built to fulfill this task. Its mathematical model works as follows 
fl2l . Satellites send to earth their location — see Figure [TJ for 
illustration. 




Fig. 1. Satellites communicate location in GPS. 

For simplicity, the location of a satellite is a bit b g {±1}- The 
satellite transmits to the earth its sequence S £ % of norm one 
multiplied by its location b. We assume, for simplicity, that the 
sequence travels through only onepath. Hence, the client receives 
the sequence R e R of the formj 

R[n] = b ■ a ■ e*&"°- n ■ S[n + t ] +W[n], (I-A.l) 

where ato E C is the complex amplitude, with |ao| < 1, G Zjv 
encodes the radial velocity of the satellite with respect to the 
client, t G 1*n encodes the distance between the satellite and 
the clienj^] and W is a random white nois^j The problem of GPS 

'We denote i = y—X. 

2 Using to we can compute 1121 the distance from the satellite to the client, 
assuming a line of sight between them. 

3 In this paper, a random white noise will always be assumed to have mean 
zero. 
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can be formulated as follows: 



Problem I-A.l (The GPS Problem): Design S G H, and an 
effective method of extracting (b, tq) from S and R satisfying 
dTXT) . 

In practice, the satellite transmits S = Si + b ■ S2, where 

51 , S2 are almost orthogonal in some appropriate sense. Then 
(a ,T ,Wo), and (b ■ a ,T ,w ) are computed using Si, and 

52, respectively, concluding with the derivation of the bit b. A 
client can compute his/her location by knowing the locations of 
at least three satellites and distances to them. The GPS problem 
is an example of channel estimation task. We would like now to 
describe the more general channel estimation problem that we are 
going to solve. 

B. Channel Estimation Problem 

We consider the following mathematical model of time- 
frequency channel estimation |16|. There exists a collection of 
users, each one holds a sequence from H known to a base 
station (receiver). The users transmit their sequences to the base 
station. Due to multipath effects — see Figure [2] for illustration — 
the sequences undergo fT31 . Ifl6l several time-frequency shifts 
as a result of reflections from various obstacles. We make the 
standard assumption of almost-orthogonality between sequences 
of different users. Hence, if a user transmits S E H, then the base 
station receives R E T~L of the form 

m 

R[n] = ^areT^".S[n + T fc ]+W[n], n € Z N , (I-B.l) 
fe=i 

where m denotes the number of paths the transmitted sequence 
traveled, E C is the complex multipath amplitude along path k, 
with Y^k=i l a fc| 2 — 1j u k € 1>n depends on the relative velocity 
along path k of the transmitter with respect to a base station, 
TjeZiv encodes the delay along path k, and W E H denotes a 
random white noise. The parameter m will be called the sparsity 
of the channel. The objective is: 

Problem I-B.l (The Channel Estimatiom Problem): Design 
S € H, and an effective method of extracting the channel 
parameters (ctfc, Tfc, Wfc), k = 1, ...,m, from S and R satisfying 



To suggest a solution to Problem I-B.l we start with a simpler 
variant. 

C. The Time-Frequency Shift (TFS) Problem 

Suppose the transmitter and the receiver sequences S,R E H 
are related by 

R[n] =e^ u °- n ■ S[n + T ] +W[n], n€Z N , (I-C.l) 

where W E H denotes a random white noise, and (to,u>q) E 
Zjv x Zjy- The pair (to,Wq) is called the time-frequency shift, 
and the vector space V = Zn x Zn is called the time-frequency 
plane. We would like to solve the following: 

Problem I-C.l (Time-Frequency Shift (TFS)): Design S E 
H, and an effective method of extracting the time-frequency shift 



(tq,Wq) from S and R satisfying (I-C.l 




Transmitter 



Fig. 2. Three paths scenario. 



Receiver 



D. The Matched Filter (MF) Algorithm 

A cla ssical solution 0, (8), 0, CU, (16), E), ED to 
Problem I-C.l is the matched filter algorithm. We define the 



following matched filter (MF) matrix of R and S: 

M(R,S)[t,uj} = (i?[rj],e^-™ -S[n + r]\, (t,lu) E V. 



A direct verification shows that for ( 
probability one, we have 



e ^i(TWo-wTo) with 



M(R, S)[t,lo] 



C -M(S,S)[t 



T ,UJ - LO \ 



N 



where NSR 



SNR' 



.e., essentially^] the inverse of the signal- 
to-noise ratio between the sequences S and W. For simplicity, 
we assume, for the rest of the paper, that the NSR is of size 
0(1). In order to extract the time-frequency shift (tq,u>q), it is 
"standard'0(see @, 0, 0, flj], ED, ED, ED) to use pseudo- 
random sequence SeHof norm one. In this case A4(S, S)[t 
t , (J - cj ] = 1 for (t,oj] 
(t,uj) ^ (r , Wo)- Hence, 



(to, w ), and of order 0(^=) if 



M(R,S)[t,u] = 



where e N = 0(-^)- 



1 + e N , if (t,w) = (t , w ); 
£at, if (r,w) ^ (t , w ), 



(I-D.l) 



Identity (I-D.l 1 — see Figure [3] for a demonstration — suggests 
the following "entry-by-entry" solution to TFS problem: Com- 
pute the matrix Ai(R, S), and choose (to,wo) for which 
A4(R, S)[tq,ojq] w 1. However, this solution of TFS problem 
is significantly expensive in terms of arithmetic complexity, i.e., 
the number of multiplication and addition operations is 0(N 3 ). 
One can do better using a "line-by-line" computation. This is due 
to the next observation. 



4 The precise relation is NSR 
logarithm. 



For example in spread-spectrum communication systems. 
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Delay 



Fig. 3. \M(R, S)\ with pseudo-random S, and (to,dq) = (50,50). 



Fig. 4. \M(R,S L )\ for a flag S L with L = {(0,^)}, and (t ,w ) = (50,50). 



I-D.l one can solve TFS 



Remark I-D.l (FFT): The restriction of the matrix M(R, S) 
to any line (not necessarily through the origin) in the time- 
frequency plane V, is a convolution that can be computed, using 
the fast Fourier transform (FFT), in (3(Alog A) operations. For 

details see Section [V] 

As a consequence of Remark 
problem in 0(N 2 log N) operations. 



E. The Fast Matched Filter (FMF) Problem 

To the best of our knowledge, the "line-by-line" computation 
is also the fastest known method fl4l . If N is large this may not 
suffice. For example in applications to GPS as in Problem 



I-A.l above, we have N > 1000. This leads to the following: 

Problem I-E.l (The Fast Matched Filter Problem): Solve 
the TFS problem in almost linear complexity. 

Note that computing one entry in A4 (R, S) already takes 
O(N) operations. 

F. The Flag Method 

In this paper we introduce the flag method to propose a solution 
to FMF problem. The idea is, first, to find a line on which the 
time-frequency shift is located, and, then, to search on the line 
to find the time-frequency shift. We associate with the A + 1 
lines Lj, j — 1, ...,N + 1, through the origin in V, a system 
of "almost orthogonal" sequences Sl- G T~L, that we call flags. 
They satisfy — see Figure [4] for illustra tion — - the following "flag 

with S = SL, , we 



property'^] For a sequence R given by (I-C.l 
have 



In addition, for Sl and R satisfying ( I-F.l| >, we have the 
following search method to solve FMF problem: 

Flag Algorithm 

Step 1. Choose a line L 1 - transversal to L. 

Step 2. Compute M(R,Sl) on L^. Find (r, uj) such that 
\M(R, Sl)[t, uj}\ s» 1, i.e., (r,w) on the shifted line 
L + (r , w ). 

Step 3. Compute M(R, Sl) on L + (tq,ujq) and find (r, ui) 
such that \M(R,S l )[t,u}\ « 2. 

The complexity of the flag algorithm — see Figure [5] for a 
demonstration — is O(NlogN), using the FFT. 



i 



L 1 



i 



Fig. 5. Diagram of flag algorithm. 



This completes our solution of Problem I-E.l 
Matched Filter Problem. 



-The Fast 



M{R,Sl s )[tM (I-F.l) 

2 + e N , if (t,u}) = (t , w ); 
l + e N , in |-| if (t,w) e L'j \ (t ,w ); 
e N , if (t, w) G V n L'j, 

where en = 0(^J=), \-\ denotes absolute value, and Lj is 
the shifted line Lj + (t ,u>q). The "almost orthogonality" of 



sequences means 

(t,w), i ^ j. 



\M(S Li ,S Lj )[r,u;}\ = 0(^=), for every 



G. Solution to the GPS and Channel Estimation Problems 
Let L C V be a line through the origin. 

Definition l-G.l (Genericity): We say that the points 
(Tfc,u;fc) € V, k = 1, ...,m, are L-generic if no two of them lie 
on a shift of L, i.e., on L + v, for some »6F. 



6 In linear algebra, a pair 
j 6 L, is called a flag. 



&0,L) consisting of a line L <Z V, and a point 
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Looking back to Problem |I-B.1| we see that, under genericity 
assumptions, the flag method provides a fast computation, in 
0(m ■ NlogN) operations, of the channel parameters of channel 
with sparsity m In particular, it calculates the GPS parameters — 
see Problem I-A.l — in 0(N log N) operations. Indeed, Identity 



(I-F.li, together with the almost orthogonality between flag 



sequences, implies that 



a k ^M{R,S L )[T k ,u k }/2, k = l, 



where R is the sequence ( |I-B.1[ ), with S = Sl, assuming that and 
(rfc,u;fc)'s are L-generic. So we can adjust the flag algorithm as 
follows: 

• Compute A4(R,Sl) on L 1 - . Find all (t,oj)'s such that 
\A4(R, Sl)[t,ui] \ is sufficiently large, i.e., find all the shifted 
lines L + (r fe ,o; fe )'s. 

• Compute A4(R,Sl) on each line L + (T k ,u) k ), and find 
(t,w) such that \A4(R, Sl)[t, oj]\ is maximal on that line, 
i.e., (t,cj) = (r fc ,cj fe ) and a k » M(R, S L )[r k , ui k ]/2. 

Figure [6] provides a visual illustration for the matched filter 
matrix in three paths scenario. 




the form 

R[n] =^2a k -e^"- n ■ S[n + T k ]+W[n], nel 



N ■ 



k=l 



where a k G C is the complex multipath amplitude along path k, 
with 2~2T=i l a fc| 2 < 1) w fc G 2jv encodes the radial velocity of 
target k with respect to the radar, r k G encodes the distance 
between target k and the radar, and W is a random white noise. 

In order to determine the location of the targets we need to 
solve [13] the following. 

Problem I-H.l (The Radar Problem): Having R and S, com- 
pute the parameters (r k ,oj k ), k = 1, ...,m. 

This is essentially the channel estimation problem. Under the 
genericity assumption, the flag method solves it in 0(m-N log N) 
operations. This completes our solution to Problem I-H.lj — The 
Radar Problem. 




Sender/ 
Receiver 



[ Object 



Delay 



distance 



Fig. 7. Radar transmits wave and recieves echo. 



Fig. 6. \M{R,S L )\, for L = {(0,w)}, and (a k ,T k ,w k ) = (i, 50fc, 50fe), 
fc = l,2,3. 



This completes our solutions of Problem I-B.l — The Channel 



Estimation Problem, and of Problem |I-A.lj — The GPS Problem. 



H. Applications to Radar and Mobile Communication 

The flag method provides a significant improvement over the 
current channel estimation algorithms in the presence of high 
velocities. The latter occurs in systems such as GPS, radar, and 



mobile communication of fast moving users. In Subsection I-A 



we described the GPS problem, and in Subsection |I-G| its effective 
solution using the flag method. It is easy to see that the flag 
method suggests a solution to the GPS problem also in the 
multipath scenario. In this section we demonstrate application of 
the flag method to radar, and mobile communication. 

1) Application to Radar: The model of radar works as 
follows ifTTl . A radar transmits — Figure [7] illustrates the case of 
one target — a sequence S € H which bounces back from m 
targets. The sequence R € T-L which is received as an echo has 



2) Application to Mobile Communication: The model of 
mobile communication works as follows 1161 . A user wants to 
deliver a bit of information b € {±1} to a base station. The 
base station assigns a sequence S & T-L to the user, and the user 
transmits to the base station the sequence b ■ S. The sequence 
R G % which is received by the base station is of the form 



R[n] = b ■ a k ■ e* 

k=l 



S[n + T k ]+W[n], neZ N , 



where m denotes the number of paths the transmitted sequence 
traveled, a k € C is the multipath amplitude along path k, with 
SfcLi \ a k\ 2 < 1, Wfc € depends on the relative velocity along 
path k of the user with respect to the base station, r k £ Ijq is the 
delay along path k, and W G H denotes a random white noise. 
The main task at the base station is the following: 

Problem I-H.2 (The Mobile Communication Problem): 
Having R and S, compute the bit b. 

In practice, first the user sends S, and the channel estimation 
is done. Then the bit b is communicated by sending b ■ S. Finally, 
knowing the channel parameters (a k ,T k ,cu k ), k — l,...,m, the 
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bit is extracted using the formula^] 

rn rn 

& • ^ ki 2 ~ (r, a * • e ^ k ' n ■ s ^ n + 



fc=i 



The main computational step is the channel estimation which is 
done by flag method in 0(m-N log N) operations. This completes 
our solution to Problem II-H.21 — The Mobile Communication 
Problem. 

/. What you can find in this paper 

• In Section |7J You can read about the flag method for effec- 
tive delay-Doppler channel estimation. In addition, concrete 
applications to GPS, radar, and mobile communication are 
discussed. 

• In Section [//]• You can find the definition and explicit formu- 
las for the Heisenberg and Weil operators. These operators 
are our basic tool in the development of the flag method, in 
general, and the flag sequences, in particular. 



In Section III You can see the design of the Heisenberg- 
Weil flag sequences, using the Heisenberg-Weil operators, 
and diagonalization techniques of commuting operators. In 
addition, the investigation of the correlation properties of the 
flag sequences is done in this section. These properties are 



formulated in Theorem III-C.l which guarantees applicabil- 



ity of the Heisenberg-Weil sequences to the flag method. 



In Section IV You can get explicit formulas for large col- 
lection of the Heisenberg-Weil flag sequences. In particular, 
these formulas enable to generate the sequences using low 
complexity algorithm. 

In Section [V| You can find the formulas that suggest fast 
computation of the matched filter matrix on any line in 
the time-frequency plane. These formulas are of crucial 
importance for the effectiveness of the flag method. 



In Section VI You can find needed proofs and justifications 



for all the claims and formulas that appear in the body of 
the paper. 
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II. The Heisenberg and Weil Operators 



The flag sequences (see Subsection I-Fi are defined, con- 
structed and analyzed using two special classes of operators that 
act on the Hilbert space of digital sequences. The first class 
consists of the Heisenberg operators and is a generalization of the 

7 It is analogous to data modulation using a delay-Doppler rake receiver in 
spread-spectrum systems [16]. 



time-shift and frequency-shift operators. The second class consists 
of the Weil operators and is a generalization of the discrete Fourier 
transform. In this section we recall the definitions and explicit 
formulas of these operators. 

A. The Heisenberg Operators 

The Heisenberg operators are the unitary transformations that 
act on the Hilbert space of digital sequences by 

7t(t, oj) : H —r H, t,oj E Zjv; 
[w{r,oj)f}[n] = e a r«-«./[n + r], 

for every f £ H, n £ In- 



B. The Weil Operators 

Consider the discrete Fourier transform 

DFT :H->H, 

N-l 

[DFT(f)M = ^E 



(II-A.1) 



71=0 



fir 



for every f e H, u E Z N . It is easy to check that DFT satisfies 
the following A 2 identities: 



DFT o 7t(t, oj) = 7r(— cj, t) o DFT, t,loeZ 



N- 



(II-B.l) 



where tt(t,uj) are the Heisenberg operators, and o denotes com- 
position of transformations. A version of the celebrated Stone-von 
Neumann (S-vN) theorem implies that up to scalar multiple the 



DFT is the unique operator that satisfies (II-B.l I. This means that 



( II-B. 1 1 is a characterization of the DFT. In [ 19 1 Weil generalized 
this method and defined many other operators that act on H. 
Consider the following collection of matrices 



SL 2 (Z 



N) 



a b 
c d 



a, 6, c, d S Zjv, and ad — be — 1 



Note that G = SL2 (Zjv) is a group |2| with respect to the 
operation of matrix multiplication. It is called the special linear 
group of order two over 7L N . Each element 



9 



e G, 



acts on the time-frequency plane V 
of coordinates 



Zjv x Zjv via the change 



(r, oj) H> g ■ (r, oj) = (ar + bo:, cr + doj). 

For g E G, let p(g) be a linear operator on H which is a solution 
of the following system of A 2 linear equations: 

s s : p{g)°^{ T ^) = 7r (5-( r , w ))°p(3), r,cj € Zjv, (II-B.2) 



Denote by Sol(S 9 ) the space of all solutions to System (II-B. 2 1. 
For example for 

'0 -1 



1 







which is called the Weyl element, we have by (II-B.l 1 that DFT g 
Sol(E w ). The general version of the S-vN theorem implies that 
dimSol(S 9 ) = 1, for every g G G. In fact there exists a special 
set of solutions. This is the content of the following result fl9l : 
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Theorem II-B.l (Weil operators): There exists a unique col- 
lection of solutions {p(g) £ Sol(E 9 ); g £ G}, which are unitary 
operators, and satisfy the homomorphism condition p(g ■ h) — 

p(g) o p(h), for every g,h £ G. 

Denote by U(H) the collection of all unitary operators on the 



Hilbert space H of digital sequences. Theorem II-B.l establishes 
the map 

p:G^U(H), (II-B.3) 

which is called the Weil representation [19|. We will call each 
p(d)i 9 S C, a Weil operator. 

1) Formulas for Weil Operators: It will be important for our 
study to have the following 0, explicit formulas for the Weil 
operators: 



Fourier. We have 



Chirp. We have 

1 

y c 1 

Scaling. We have 

a 
a- 1 



DFT(f) [r 



K-2-W) . f[r 



f 



N 



f[a 



(II-B.4) 



(II-B.5) 



(II-B.6) 



for every / g H, 0^fl,c,ne Zjv, where N^J is the Legendre 
symbol which is e qual to 1 if a is a square modulo N, and —1 
otherwise, and in (II-B.5 1 we denote 2 _1 = =^±1 the inverse of 
2 modulo N. 

The group G admits [3] the Bruhat decomposition 

G = U A U UwUA, 

where U C G denotes the unipotent subgroup 

'1 



U 



c e Z 



and A C G denotes the diagonal subgroup 



A = 







a 



0^a£ 



(II-B.7) 



This means that every element g £ G can be written in the form 

g = u ■ s or g = u ■ w • u ■ s 

where u, u' , u" £ U, s, s' € A, and w is the Weyl element. Hence, 
because p is homomorphism, i.e., p(g ■ h) — p(g) o p(h) for every 
g, ZieG.we deduce that formulas ( |II-B.4| ), ( |II-B.5| l, and ( |II-B.6| l, 
extend to describe all the Weil operators. 

III. Sequence Design: Heisenberg-Weil Flags 

The flag sequences, that play the main role in the flag method, 
are of a special type. Each of them is a sum of a pseudorandom 
sequence and a structural sequence. The first has the MF matrix 
which is almost delta function at the origin, and the MF matrix of 
the second is supported on a line. The design of these sequences 
is done using group representation theory. The pseudorandom 
sequences are designed [8|, [9], [18| using the Weil representation 



operators (II-B.3 I, and will be called Weil (spike) sequence^ The 
structural sequences are designed ifTUll . ifTTIl using the Heisenberg 
representation operators ( II-A. 1| >, and will be called Heisenberg 
(line) sequences. We call the collection of all flag sequences, the 
Heisenberg-Weil flag system. In this section we study construc- 
tions, and properties of these sequences. 



A. The Heisenberg (Lines) System 

The operators ( II-A. 1[ > obey the Heisenberg commutation 
relations 



w(t,u>) o7t(t' ,J) = eT^™'-^') ■7r(T',u/)°'r(^)- 

The expression tlu' — lot' vanishes if (r, tS), (t',u>') are on the 
same line through the origin. Hence, for a given line L C V = 
Zjv x Zjv we have a commutative collection of unitary operators 



tt(£) :H^H, £eL. 



(III-A.l) 



Explicit version of simultaneous diagonalization theorem from 
linear algebra implies iflOl . ifTTl the existence of a natural or- 
thonormal basis Bl for % consisting of common eigensequences 
for all the operators ( III-A. 1 [ > 

where ip runs over characters of L, i.e., functions ip : L — > C* = 
C - 0, with ip{i + I') = ip[f)<ip{l'), for every £,£' e L. The 
system of all such bases Bl , where L runs over all lines through 
the origin in V, will be called the Heisenberg (lines) system. We 
use the following result |10|, OH: 

Theorem III-A.l: The Heisenberg system satisfies the follow- 
ing properties: 

1) Line. For every line L C V, and every /l £ Bl, we have 



M{f L ,f L )[rM 



1, if (r,«) = (0,0); 

1, in |-| if (t,uj) £ L \ (0,0); 

0, if (t,u>) i L. 



2) Almost-orthogonality. For every two lines L\ ^ L2 C V, 
and every f Ll £ B Ll , Jl 2 £ Bl 2 , we have 

1 

for every (t, u) £ V. 



\M(f Ll ,f L2 )[T,u;] 



Figure [8] demonstrates Property 1 of Theorem III-A. 1 for the 
diagonal line. 

B. The Weil (Spikes) System 

The group G = SL2 (Zjv) is non-commutative, but contains a 
special class of maximal commutative subgroups called torj^]|8|, 
1 9 1, |3|. Each torus T C G acts via the Weil operators 



p(g):H^H, g £ T. 



(HI-B.l) 



8 For the purpose of the Flag method, other pseudorandom signals may work. 
There are order of N 2 tori in SL 2 (Z N ). 
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Delay 




Delay 



Fig. 8. \M(f L , f L )\ for L = {(t,t))}. 



Kg. 9. X[ Vt , Vt ] forT= {I Q a _i 1 ;0^aGZjv} 



a 



This is a commutative collection of diagonalizable operators, and 
it admits [8 1, [9] a natural orthonormal basis Bt for W, consisting 



of common eigensequences for all the operators (III-B.l 



B T = {ip Tx }; 



(HI-B.2) 



where x mns over characters of T, i.e., functions \ : T 
with X (g ■ 9 1 ) = X(9)x{g'), for every g l9 ' G T. 



Remark III-B.l: There is a small abuse of notation in IIII-B.21 
There are two types of tori in G, split tori and non-split tori 
|H), 0. In each case, the torus T admits a unique non-trivial 
character \ q of T — called the quadratic character — which takes 
the values Xgid) G {^l}; 5 € T. The dimension of the space 
n x of sequences <p T , which satisfy p(g)tp T = X„{g)^ is 

y Aq Xq ^ X. n 

equal to 2 or 0, if T is a split or non-split torus, respectively [8|, 
0. 



Let us denote by 



the set of sequences in Bt, which are not associated with the 
quadratic character. The system of all such sets St, where T 
runs over all tori in G, will be called the Weil (spikes) system. 
We use the following result Qj), 0: 

Theorem III-B.2: The Weil system satisfies the following 
properties: 

1) Spike. For every torus T C G, and every ip T G St, we 
have 



M(ip T ,tp T )[r,ui\ 



1, if (t,w) = (0,0); 



< 



/AT 



in 



if (t.wJ^O.O). 



2) Almost-orthogonality. For every two tori Ti, T 2 C G, and 
every yj T G 5^ , <Pt 2 G St 2 , with <^ Ti ^ (p T2 , we have 



|-A^(9? Tl ,(/3 r2 )[T,w]| < 
for every (r, uS) G V. 



/AT 



/JV ' 



if Ti ^T 2 ; 
if T, =T 2 , 



Figure [9] illustrates Property 1 of Theorem III-B.2 applied to 
the commutative subgroup of diagonal matrices in G. 



C. The Heisenberg-Weil System 

We define the Heisenberg-Weil system of sequences. This is 
the collection of sequences in H, which are of the form Sl = 
II + <Pt> where Jl and <p T are Heisenberg and Weil sequences, 
respectively. The main technical result of this paper is: 

Theorem III-C.l: The Heisenberg-Weil system satisfies the 
properties 

1) Flag. For every line L C V, torus T C G, and every flag 
Sl = II + <Pt> with f L e B L , ip T G S T , we have 



M(S l ,S l )[t,lo] 



2 + e N , if (r,w) = (0,0); 

1 + £jv, in |-| if (t,w) G L \ (0, 

£at, if (r,w) eF\L, 



where I ejv I < 



hf> and 1^1 ^ > 



2) Almost-orthogonality. For every two lines L\ ^ T 2 C V, 
tori Ti, T 2 C G, and every two flags = f Lj + ip T ., 
with f L . G S Lj , p T G S Tj , j = 1, 2, </? Tl 7^ <^ T2 , we have 



\M(S Li ,S L2 )[tM\<\ ^ 



I N 



if Ti ^T 2 ; 
if Ti = T 2 , 



for every (t, cj) G V. 



For a proof of Theorem III-C.l see Subsection VI- A 



Figure 10 demonstrates Property 1 of Theorem III-C.l applied 



to the diagonal line L and a torus T as in Figure 9 



Remark III-C.2: As a consequence of Theorem III-C.l we ob- 



tain families of N+l almost-orthogonal flag sequences which can 
be used for solving the TFS problem in 0(N\ogN) operations, 
and GPS, channel estimation, radar, and mobile communication 
problems in 0(m ■ NlogN) operations for channel of sparsity 
m. 

This completes our design of the Heisenberg-Weil flag se- 
quences. 
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Fig. 10. \M[S L ,S L ]\ for Heisenberg-Weil flag with L = {{r, r)}. 

IV. Formulas for Heisenberg-Weil Sequences 

In order to implement the flag method it is important to have 
explicit formulas for the Heisenberg and Weil sequences, which 
in particular enable one to generate them with a low complexity 
procedure. In this section we supply such effective description 
for all Heisenberg sequences, and for Weil sequences associated 
with split tori 

A. Formulas for Heisenberg Sequences 

First we parametrize the lines in the time-frequency plane, and 
then we provide explicit formulas for the orthonormal bases of 
sequences associated with the lines. 

1 ) Parametrization of Lines: The N + 1 lines in the time- 
frequency plane V = In x Zjv can be described in terms of their 
slopes. We have 

• Lines with finite slope. These are the lines of the form 
L c =span{(l. c)}, c £ Zjv. 

• Line with infinite slope. This is the line L^ =span{(0, 1)}. 

2) Formulas: Using the above parametrization, we obtain 

• Formulas for Heisenberg sequences associated with lines of 
finite slope. For c £ Z^r we have the orthonormal basis 



Bl c = {/ c ,bN = -j=e 
V A 



of Heisenberg sequences associated with the line L c . 

Formulas for Heisenberg sequences associated with the line 
of infinite slope. We have the orthonormal basis 

B Lao = {S b ; b £ Zat}, 

of Heisenberg sequences associated with the line L^ , where 
the Jb's denote the Dirac delta functions. 



1) Parametrization of Split Tori: A commutative subgroup 
T C G is called split torus Q if for some g £ G it is of the 



form T — T g , with 



T g =g-A-g-\ 



where A C G is the subgroup of all diagonal matrices, also called 
the standard torus, i.e., 

A = {(l -i); 0^aeZ N 

We denote by T = {T g ; g € G} the set of all split tori in G. It is 
not hard to verify that the number of elements in T is iV (-^+ 1 ) , 
A direct computation shows that the collection of all T 9 's with 



(IV-B.l) 



exhausts the set T. Moreover, in (IV-B.l I the torus T g can be 
written also as T g r, for g ^ g' , only if b ^ and 
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1 



1 + be 







2) Formulas: In order to provide the explicit formulas we 
need to develop some basic facts and notations from the theory 
of multiplicative characters [2 |. Consider the group Z* N of all non- 
zero elements in Z^ , with multiplication modulo N. A basic fact 
about this group is that it is cyclic, i.e., there exists an element 



r e Z* N such that 



Z 



JV 



{l,r,r 2 



„N-2 



}• 



A function \ '■ ~ * ^* ^ s called multiplicative character if 
x(x ■ y) — x( x ) ' x(y) f° r every x, y <E Z* N . A way to write 
formulas for such functions is the following. Choose ( e C which 



satisfies £ 



7V-1 



l,i.e., Ce^jv-i ={e^ k ; k = 0,...,N -2}, 



and define a multiplicative character by 



X C (r d ) = C 



0,1, ...,7V -2. 



Running over all the — 1 possible such £'s, we obtain all the 
multiplicative characters of Z* N . We are ready to write, in terms 



of the parametrization (IV-B.l I, the concrete eigensequences 



associated with each of the tori (see Subsection |HI-B| >. We obtain 

• Formulas for Weil sequences associated with the diagonal 
torus. For the diagonal torus A we have 



where (p x € H is the sequence defined by 



if n = 0. 



(IV-B.2) 



B. Formulas for the Weil Sequences 

We describe explicit formulas for the Weil sequences associ- 
ated with split tori 0, QH], 0. First we parametrize the split tori 
in G = SX 2 (Zjv), and then we write the explicit expressions for 
the orthonormal bases of sequences associated with these tori. 



Formulas for Weil sequences associated with the torus T Uc , 
for unipotent u c £ G. For the torus T Uo associated with the 
unipotent element 



1 

c 1 



c £ Z 
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we have 



where <p x ^c £ % is the sequence defined by 



e N 



(IV-B.3) 



for every n £ Zjv, and (p x is the sequence given by 

Hrv-B.2) . 

Formulas for Weil sequences associated with the torus T gi 
for non-unipotent g £ G. For the torus T g associated with 
the element 



9 



1 b 

c 1 + bc 



b, c £ Zj\r, b 7^ 0, 



we have 



where ip x s is the sequence defined by 



= C& 



e « 



i(2 



-1 1 + bc 



E 



e w 



e N 



i(-2 _:L i>a; 2 ) 



(IV-B.4) 



ip X( [bw] 



for every n G Zjv, and <^ x ^ the sequence given by (|IV-B.2|, 

C b = i~" 



(f) 



with the Legendre symbol. 



The validity of Formula ( |IV-B.2 1 is immediate from Identity 
( |II-B.6[ ). For a verification of Formulas ( |IV-B.3| > and ( |IV-B.4| >, see 
Subsection IVFBl 

Remark IV-B.l (Complexity of Heisenberg-Weil sequences): 
For concrete applications it is important to have low arithmetic 
complexity algorithm generating the sequences. Note that the 



sequence (IV-B.4i can be computed in 0(N log N) operations 
using FFT. We conclude that the all Heisenberg sequences, and 
all Weil sequences associated with split tori, and in particular 
the associated flag sequence, can be computed in at most 
0(N log N) operations. 



1) Formula on lines with finite slope and their shifts. For 
c £ Zjv consider the line L c = {t ■ (1, c) ; r £ Zjv}, and 
for a fixed w £ Zjv the shifted line L' c = L c + (0,w). On 
L' c we have 

M(<p,4>)[r-(l,c) + (Q,u)] (V-.2) 

= [mcxp(2- 1 cn 2 +wn) 1 ^- * rn oxp (_ 2 - 1 cn 2 )^] [t], 

where [m exp(2 -i cn 2 +am) ip_] [n] = e 2 ^ ( 2 ~ lc " 2 + w ™) x 

<£>_[n], n £ Zjv, and similar definition for the second 
expression, with 4> the complex conjugate of the sequence 
0. 

2) Formula on the line with infinite slope and its shifts. 
Consider the line = {uj ■ (0, 1) ; uj £ Zjv}, and for 
a fixed r £ Zjv the shifted line L'^ = L x + (r, 0). On 
L'qq we have 

J M(v>,^)[w-(0 J l) + (r J Q)]=£>fT(v>-^)H, (V-.3) 
for every uj £ Zjv- 



The validity of Formula ( V-.3 1 is immediate from the definition 
of the matched filter. For a verification of Formula ( |V-.2| i see 
Subsection IVFCl 



CL 
O 



Delay 



Fig. 11. M(<p,4>)[t,t] = m cxp(2 -i n2) ^_ * m exp( _ 2 _i n 2)<£ M on L\. 



V. Computing the Matched Filter on a Line 

Implementing the flag method, we need to compute in 
0(N log N) operations the restriction of the MF matrix to any 
line in the time-frequency plane (see Remark I-D.li. In this 
section we provide algorithm that fulfills this task. The upshot 
is — see Figure 11 for illustration of the case of the diagonal 
line — that the restriction of the MF matrix to a line is a certain 
convolution that can be computed fast using FFT. Denote by 
M(<p, 4>)[t,uj] — (ip, 7r(r,w)0) the matched filter associated 
with sequences <p, cj> £ H, and by <p * <fr G H their convolution 



(ip * 4>) [t] = <P-[n] ■ fH. 



(V-.l) 



where y>_[ri] = (f[— n], and <j> T [n] = 4>[t + n], for every r,n £ 



UN- 



We consider two cases: 



VI. Proofs 



A. Proof of Theorem III-C.l 
1 ) Flag Property 

M(S L ,S L ) 



Let St = ft + <Pt- We have 

= M(f L J L ) + M(f L ,<p T ) 
+M(p T , ft) + M(ip T ,cp T ) 

We will show that 

\M(<p. T ,ft)[T,Lj}\< 



N 



t ,U) G Z 



AT- 



(VI-A.l) 



Notin g that M (f L ,ip T )[T,uj] = M(tp T , ft)[~r, -u] we get 
from ( VI-A.l \ also the same bound for M-ifttfr)- Having 
this, using Theorems |III- A. 1 1 and |III-B~2 we can deduce the Flag 
Property. So assume Tr(£)ft — i>{t)ft for t. £ L. We proceed in 
several steps: 
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Step 1. For every v £ V we have 7r(i>)_/x £ Bl- Indeed, for 
£ £ L we have 

TT(£)[7r(v)f L ] = e 2 ^- n ^ir{vy{£)J L 
= e^- n V> v »iP(£)[ir(v)f L }, 

where £1 : V x V — > Zjv is the symplectic form 
f2[(T,w), (t',uj')] = tuj'—ojt'. Namely, 7t(u)/l is eigensequence 
for w(£) with character ip v (£) = e T(- n M)^(4 
By step 1, it is enough to bound the inner product 



a) Proof of Lemma |V/-A.7| : For I £ L we have 

= mp(s)h 



where the second equality is by Identity ( |II-B.2 >. This completes 



the proof of Lemma VI-A. 1 



\(<p T ,h)\ < 



N 



(VI-A.2) 



Step 2. The bound ( VI-A.l \ holds for L^. Indeed, then = 
§b for some b £ Zjy, hence 

\{<PT>h OB )\=\<p T [b]\< SU P \<PtM\- 

In (9) it was shown that for every Weil sequence <p T we have 

2 

sup \v T [n}\ < — 



Step 3. The bound { VI-A.l) holds for every line L. We will 
use two lemmas. First, let L, M C V be two lines, and g £ G 
such that gL = {g ■ £ ; £ £ L} = M. For a character ip : L — > C*, 
define the character ip 9 : M — > C*, by ip 9 (g-£) = ip(£), for every 
£ G L. We have 

Lemma VI-A.l: Suppose /x, is a ?/)-eigensequence for L, i.e., 
7I '(^)/l = ipWfL, for every £ £ L. Then the sequence = 
p(g)fL is V- ,9_e ig ense q uence f° r -W- 



i>) Proof of Lemma VI-A.2\ : For h E T we have 



P{g ■ h ■ g~ 1 )<p Tg = p(g ■ h- g- 1 )p(g)ip T 

= p(g)p(h)<p T 
= x{h)p(g)fT 
= x g (g ■ h ■ g _1 )^ T3 , 

where the second equality is because p is homomorphism (see 



Theorem II-B.li. This completes our proof of Lemma VI-A.2 



and of the Flag Property. 

2) Almost Orthogonality : Let 5^ = /j, + (^ Tj , j = 1,2, 
as in the assumptions. We have 

M(S Ll ,S L2 ) = M(f Ll J L2 )+M(.f Ll , VT2 ) 

+M (f Tl Jl 2 )+M(<p Ti , (p T2 ) . 



The result now follows from Theorem |III-A.1| Theorem |III-B.2| 
and the bound ( |VI-A.1[ ). This completes our proof of the Almost 
Orthogonality Property, and of Theorem III-C.l 



For a proof of Lemma VI-A.l see Subsection VI-Ala 



For the second lemma, consider a torus TcG, and an element 
g E G. Then we can define a new torus T g — gTg^ 1 = {g-h-g^ 1 ; 
h £ T}. For a character \ : T — > C*, we can associate a character 
:*, by x 9 (g ■ h ■ g- 1 ) = x(h), for every h £ T. We 



X 9 --T g 
have 



Lemma VI-A.2: Suppose ip T is a x-eigensequence for T, i.e., 
p(h)(p T = x{h)(p T , for every h £ T. Then the sequence ip T ^ = 
p(g)(p T is x 9 -eigensequence for T g . 



For a proof of Lemma VI-A.2 see Subsection VI-Alb 



Now we can verify Step 3. Indeed, given a line L C V, there 
exists g £ G such that g ■ L — L^. In particular, by Lemma 
VI-A.l we get that /l x = p{g)fh is up to a unitary scalar in 
Bl^ ■ In addition, by Lemma VI-A.2 we know that ip T — p(g)ifix 



is up to a unitary scalar in Br g - Finally, we have 

(<PT>h) = (p(g)<PT>p(g)fL) 

where the first equality is by the unitarity of p(g). Hence, by Step 
2, we get the desired bound also in this case. 



B. Verification of Formulas ( IV-B.3) , and ^V-B.4) 



The idea is to use the fact that for g £ G, the explicit Weil 
operator p(g) maps the explicit set Sa to the set St , T g = 
gAg~ Y . In details, for a character x '■ A ~ ► C* and an element g £ 
G define the character x 9 : T g -> C*, by x 9 {g • h ■ g^ 1 ) = x{h), 
for every h £ A. Using Lemma 



VI-A.2 



we deduce that if tp £ 
Sa is eigensequence of A with character \, then (p x <, — p(g)ip x £ 
St is eigensequence of T g with character x 9 ■ Specializing to the 
character X = Xqi ~ 1^(6 Mjv-I; °f A, and the associated 
sequence ip x 
the formulas. 



£ Sa given by (IV-B.2i, we can proceed to verify 



1) Verification of Formula ( IV-B.3) : For the unipotent ele- 
ment 

'1 X 



c 1 



c£Z 



N- 



we have 



p(u c )ip 



e N 



where the second equality is by Formula ( II-B.5[ >. This completes 
our verification of Formula ( |IV-B.3| l. 
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2) Verification of Formula (IV-B.4<: For the element 

'\ b 
,c l + bc. 



b,ceZ N , b^O, 



its Bruhat decomposition is 



1 b 

c l + bc 



1 



1 

-1 



1 

b 1 



b- 1 N 
b y 
(VI-B.l) 



This implies that for n e Zjv we have 

ip x 9 [n] = 



e « 



(-2- 



= ™ 2 ) 



iV 



E 



> e 



with C/, = 



/, 2 where in the second equality we use 



identity ( VI-B.1[), the fact that p is homomorphism, and the For- [7 



mulas (II-B.4 1, (|II-B.5 1, (II-B.61. This completes our verification 
of Formula dIV-B.4). 



C. Verification of Formula ( V-.2) 



We verify Formula (V-.2i for the matched filter M.(<p, 4>), 
(p,4> G H, restricted to a line with finite slope. We proceed in 
two steps: 

Step 1. The formula holds for the line Lq and its shifts. Indeed, 
for a fixed u € Z^r we compute the matched filter on L' = 

L Q + (0,w) = {(r,u); r € Z/v}. We get 



M((p,<p)[r,oj] 



((p, tt(t,u)4>) 



e -^r w - n - (pin], (f)[n + T 
= |_ m exp(am)^- * </>] M, 



where the fourth equality is by the definition ( | V-. 1 1 > of *, and the 
definition (V-.2i of m exp (.). 



Step 2. The formula holds for the lines L c , c € Z^v, and their 
shifts. Indeed, the element 



1 

-c 1 



€ G, 



satisfies 



U_e-(1,C) = (1,0), 

u_ c -(0,w) = (Q,w). 



(VI-C.l) 



For a fixed to € Zjv we compute the matched filter on L^, = 

L c + (0, w) = {r • (1, c) + (0, w); r € Z^}. We get 

X(^,0)[r-(l,c) + (O, W )] 
= ( v> 7r[T.(l,c) + (O,w)]0) 
= (p(u-c)tp, p(u- c )tt[t ■ (l,c) + (0,w)]</)) 
= (p(U- c )<P, k(t,lj)p(u_ c )<P) 
= M(p(u-c)<p,p(u- c )<i>)[T,uj] 

= [ m c^p(2~ 1 cn 2 +un) l P- * m exp(-2- 1 cn 2 ) </>] H, 



where, the second equality is by the unitarity of p, the third 
equality is by Identities (II-B.2I, (VI-C.li, and the last equality 



[6 



is by Formula (II-B.5i and Step 1 above. 

This confirms Step 2, and completes our verification of For- 
mula {VT2]). 
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